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ABSTRACT 

X-ray spectra from stellar coronae are reprocessed by the underlying photo- 
sphere through scattering and photoionization events. While reprocessed X-ray 
spectra reaching a distant observer are at a flux level of only a few percent of 
that of the corona itself, characteristic lines formed by inner shell photoioniza- 
tion of some abundant elements can be significantly stronger. The emergent 
photospheric spectra are sensitive to the distance and location of the fluoresc- 
ing radiation and can provide diagnostics of coronal geometry and abundance. 
Here we present Monte Carlo simulations of the photospheric Kai,a2 doublet 
arising from quasi-neutral Fe irradiated by a coronal X-ray source. Fluorescent 
line strengths have been computed as a function of the height of the radiation 
source, the temperature of the ionising X-ray spectrum, and the viewing angle. 
We also illustrate how the fluorescence efficiencies scale with the photospheric 
metallicity and the Fe abundance. Based on the results we make three com- 
ments: (1) fluorescent Fe lines seen from pre-main sequence stars mostly suggest 
flared disk geometries and/or super-solar disk Fe abundances; (2) the extreme 
^ 1400 mA line observed from a flare on V 1486 Ori can be explained entirely 
by X-ray fluorescence if the flare itself were partially eclipsed by the limb of the 
star; and (3) the fluorescent Fe line detected by Swift during a large flare on 
II Peg is consistent with X-ray excitation and does not require a coUisional ioni- 
sation contribution. There is no convincing evidence supporting the energetically 
challenging explanation of electron impact excitation for observed stellar Fe Ka 
lines. 

Subject headings: Sun: X-rays — methods: numerical — stars: coronae — stars: 
flare — stars: pre-main sequence — X-rays: stars 
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Introduction 



It is well-established from surveys of the sky at EUV and X-ray wavelengths that all 
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outer atmospheres akin to that of the Sun (e.g. 
much observational and theoretical effort has been devoted to understanding solar coro- 
nal spectra and, in more recent years, toward understanding stellar coronal emission and 
spectra, comparatively little attention has been devoted to the reprocessing and line flu- 
orescence resulting from this coronal emission by the underlying solar and stellar photo- 
spheres. In contrast, considerable effort has been spent on the study of X-ray reprocessing 
by "cold" gas in much more complex systems with more prominent fluorescent features 
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photospheric fluorescence by coronal irradiation are the same as those discussed in these 
works; the main difference here is in the specific geometry of the X-ray source above a 
quasi-neutral sphere and of the extended, shell-like nature of the coronal source above the 
photosphere during quiescent conditions. 

X-rays emitted from a hot (T > 10^ K) corona incident on the underlying photosphere 
can undergo either Compton scattering or photoabsorption events through the ionization of 
atoms or weakly ionized species. Through scattering events, photons can be reflected back 
in a direction towards the stellar surface where they have a finite chance of escape. Compton 
scattering redistributes the spectrum to lower energies by ~ E'^ jmf.(? per collision, where E 
is the photon energy and me the electron rest mass. The spectrum reflected from a stellar 
surface by scattering events is then shifted and broadened towards lower energies. Photoion- 
ization events involving X-ray photons directed toward the photosphere are predominantly 
inner shell interactions with astrophysically abundant elements, the outer and valance cross- 
sections being very small at these energies. Observable fluorescent lines can then arise as 
a result of the finite escape probabilities of photons emitted in outward directions by hole 
transitions in these atoms photoio nized in their in ner s hells . Thes e processes have been de- 
scribed in the solar context by, e.g., Tomblin] ( 1972 ) and iBai (1979) (B79), and more recently 
for arbitrarily photoionized slabs by iKallman et al. (2004). 



The strongest of the fluorescent lines for a plasma of approximately cosmic compo- 
sition is the 2s- Ip 6.4 keV Fe Ka doublet occurring following eject ion of a Ig electron . 
It has been observed in solar spectra on numerous occasions (e.g. iNeupert et al.l 119671 ; 
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Doschek et al.lll97ll : iFeldman et al.lllQSd : iTanaka et al.lll984l : iParmar et al.lll984j : IZarro et al. 
19921) ■ The mechanism o f fluorescence by the thermal X-ray coronal continuum was sug gested 



by iNeupert et al.l (Il967f). and was firnaly es tablished on more theoretical grounds by iBasko 
f ll978l . Il979h and B79. IParmar et al.l Jl984h provided convincing observational confirmation 
based on flare spectra obtained by the Solar Maximum Mission, though it has also been 
noted that contributio ns from non-thermal electron irnpact might also be present d uring 
hard X-ray bursts (e.g. lEmslie. Phillips, fc Dennislll986l : IZarro. Dennis, fc Slateiill992l ). 



B79 pointed out that, for a given source spectrum, the observed flux of Fe Ka photons 
from the photosphere depends on essentially three parameters: the photospheric iron abun- 
dance; the height of the emitting source: and the heliocentric angle between the emitting 
source and observer. iPhillips et al.l (119941 ) used Fe K/3 observations to probe the difference 
between the photospheric and coronal iron abundance for flares observed by the Yohkoh 
satellite. More importantly for the stellar case, the spatial aspects of photospheric fluo- 
rescent line formation suggest its application to understanding the spatial distribution of 
coronal structures and flares on stars of different spectral type and activity level to the Sun 
( iDrake et al.lll999l ). Inde ed, Fe K fluoresce nce has recently been detected duri ng flares on 



the a ctive binary II Peg (jOsten et al.l 120071 ) and on the single giant HR 9024 (ITesta et al. 



20071 ). The Fe K line has also been seen in a growing sample of pre-main sequence (PMS ) 



stars (llmanishi et al.ll200ll : iFavata et al.ll2005l : iTsujimoto et al.ll2005l : iGiardino et al 



2003) 



in which the line is thought to originate predominantly from the irradiated protoplanetary 
disk rather than the photosphere. 



Since the work of iBail (119791 ). there have been no concerted efforts to extend models 
of photospheric fluorescence for coronal excitation sources with other characteristics. Fluo- 
rescent lines other than Fe Ka have also not yet, to our knowledge, been studied by other 
workers in any detail in this context. A reasonably strong feature ob served in solar spectra 
near 17.62 A had be en identifled with Fe La photospheric fluorescence ( McKenzie et al. 198ol: 



Phillips et al.lll982l) . but calculations of the expected line strength was shown by lDrake et al. 
( 119991 ) to be much too weak to explain the feature, and these authors instead identifled the 
line with a transition in Fe XVIII arising from conflgur ational mixing and b oth seen in Elec- 
tron Beam Ion Trap spectra and predicted by theory (ICornille et al.lll992l ). However, given 
the potential diagnostic value of photospheric fluorescence, other lines, such as O Ka, are 
possibly observable with very high quality observations and warrant further study. 

The capabilities of current X-ray missions such as Chandra, XMM-Newton, Swift and 
Suzaku to detect fluorescent lines further motivates a re-examination of the photospheric 
fluorescence problem in the context of stellar coronae, photospheres and protoplanetary 
disks. We restrict the study in hand to Fe K fluorescence from stellar photospheres and 
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defer detailed discussions of protoplanetary disks and fluorescence from other elements to 
future work. 



2. Calculations of Photospheric Fluorescent Spectra 
2.1. Photospheric Penetration of Coronal X-rays 



In order to establish the atmospheric conditions under which fluorescence by, and Comp- 
ton scattering of, coronal X-rays take place, it is necessary to know where in a stellar atmo- 
sphere incident X-rays are absorbed. Since the absorption cross-section for a cool plasma 
with solar, or near-solar, composition is a strong function of photon energy, the absorption of 
the coronal X-ray spectrum will take place at different levels in the atmosphere. To determine 
the characteristic height of absorption as a function of photon energy we have computed the 
optical depth for X-rays as a func tion of altitude in the solar atmosphere. We adopted the 
Vernazza. Avrett. fc Loesen (jl98ll ) Model C (VALC) describing the solar atmo spheric struc 



tureas a fu nction of altitude, and use d the photoionization cross-sections from IVerner et al. 



(I1993I ) and IVerner fc YakovlevI (Il995l ) to compute the optical de pth. The total absorption 
cross-section for a plasma with the solar chemical composition of iGrevesse fc Sauvall (119981 ) 
(GS) in the energy range of interest for this study is illustrated in Figured] for the case of 
both a neutral gas and a plasma in which all species are once-ionized (this latter situation 
does not correspond to any physically valid thermal equilibrium condition, but is simply for 
comparative purposes). 

The height in the VALC model corresponding to optical depth of unity for X-rays is 
illustrated in Figure [21 Also illustrated is the atmospheric temperature as a function of 
height in this figure for comparison purposes. At the energy corresponding to the absorption 
edge of C, optical depth of unity is reached in the lower chromosphere, at temperatures of 
~ 6000 K and about 1000 km above the point where the continuum optical depth at 5000 A, 
Tsooo; is unity (defined to be at a height of km in this model). Towards higher energies, the 
decreasing photoabsorption cross-section leads to increasingly deeper penetration of coronal 
X-rays. At the Fe K photoionization threshold energy, the typical height of absorption is 
well into the photosphere and only 100 km or so above Tsooo = 1- At this point in the solar 
atmosphere, most Fe is in the form of Fe"*" and the temperature is ~ 5400 K. 

Number densities in this region of the VALC model range from n> 10^^ down to 
10^^ cm~'^ and decrease roughly ex ponentially with an scale height ~ 130 km. Thus the 
ionization parameter, U ~ F/n (e.g. [Davidson fc Netzerlll979l ). where F is the ionizing pho- 
ton flux per unit area, is small for even the strongest flares and X-ray photoionization of 
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the photospheric gas is inconsequential. This is in contrast to models of, e.g., X-ray irradi- 
ated accretion disks where the ionization state of the reflecting medium is governed by the 
incident hard photon flux and must be computed self-consistently with the emergent X-ray 



spectrum (e.g. iKallman et al.l 120041 ). 



The range of heights in the photosphere where the X-ray photons are absorbed is very 
small compared with the stellar radius (of order ~ 100-1000 km), and the variation in 
temperature and ionization is also small over this altitude range (Figure [2]). We therefore 
approximate the photosphere by a geometrically thin but optically very thick shell of homo- 
geneous density. We set the density to an arbitrarily high value such as to obtain the same 
behaviour as for a plane-parallel model where the photosphere is approximated by a cold 
semi-infinite slab with constant density. The appropriate atmospheric reference frame for 
describing the absorption, scattering and fluorescent processes is the mass or particle column 
depth and the calculations are independent of the absolute value of the density. Provided 
plane-parallel geometry remains a good approximation, these calculations should therefore 
be valid for stars with different surface gravities and atmospheric pressures. 

By assuming the photosphere to be cold, we make no distinction between fluorescence 
processes in species of different charge states. In typical late-type stellar photospheres the 
dominant species are neutral and once-ionized, for which the splitting betwee n the respec- 



tive fluoresceiice liri es is very small as are differences in fluorescent yields ( iHousd Il969 



Kallman et al.l l2004j ) . We also note that, as illustrated in Figure (H in the energy range 
of interest here there is essentially no difference in the photoabsorption cross-sections for 
neutral or once-ionized species except at low energies where the difference is dominated by 
He and H. In reality, it is only the metals with fairly low first ionization potentials that are 
significantly ionized in the upper photosphere and lower chromosphere while He and H are 
almost entirely neutral. 



2.2. Monte Carlo (MC) Calculations of Fluorescent Processes 



We use a modified version of t he MC photoionization and radiative transfer code mo- 
cassin (lErcolano et al.l l2003l . 120051 ). This code uses a stochastic approach to the transfer of 
radiation which allows the construction of models of arbitrary geometry where all compo- 
nents of the radiation field are treated self-consistently. The basic idea is to employ a discrete 
description of the radiation field, whereb y the simulation quan ta are monochromatic pack- 
ets of radiation of constant energy, Eq, (lAbbott &: Lucyl Il985l ). The individual processes 
of scattering, absorption and re-emission of radiation can then be simulated by sampling 
probability density functions based on the medium opacities and emissivities, thus allowing 
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the determination of the trajectories of packets as they leave the illuminating source(s) and 
diffuse outward. These trajectories constitute the MC observable of our experiment which 
can be rela ted to the mean intensity of the radiation field via the following MC estimator 
flLucvlll999h : 

4vrJ.rf.= ^ij^/, (1) 



du 



where V is the volume of the current grid cell and the summation is over all the fragments 
of trajectory, /, in V, for packets with frequencies in the {u, v + dv) interval. A t is the 
duration of the Monte Carlo experiment. The significance of this estimator may be most 
easily understood by considering that the energy density of the radiation field in (z/, v + 
dv) is An Jydv/c. At each given instant a packet contributes energy Eq to the volume 
element containing it. The time- averaged energy content of a volume element is given by the 
summation of the energy contribution of all packets crossing this volume element in the time 
interval At. Equation 1 follows from this argument, since each packet contributes 6^51/ At 
to the time averaged energy content of a volume element and 5t = l/c, where I is the segment 
of the trajectory of a packet contained in the volume element. 

mocass i n's at omic database is regularly updat ed and currently uses opacity data from 
Verner et al.l (119931 ) and IVerner fc YakovlevI (119951 ). energy levels, collision strengths and 
transition probabilities from Version 5 of the CHIANTI database (ILandi et al.l |2006| . and 
references therein) and the improved hydrogen and helium free-bound continuous ern ission 
data recently published by lErcolano &: Storeyl (120061 ). The IVerner &: YakovlevI (119951 ) pho- 
toabsorption cross-sections do no not include the resonance structure that can be prominent 
near ionisation t hresholds. While this structure can be important for some astrophysical 
applications (e.g. iKallman et al.l 120041 ). it is not significant for the fluorescent efficiencies 



discussed here that depend primarily on broad-band spectral properties and neutral or once- 
ionised gas. 

In the case of inner shell FeKa production, the fluorescent luminosity, L(FeKa), due 
to the absorption of an incident energy packet of wavelength shorter than the Fe K edge 
{hv > 7.11 keV) is calculated on the fly. The event of absorption of a high energy packet 
is immediately followed by the re-emission of a number, n of Fe Ka packets from the same 
event location. For an incident packet of frequency z/, carrying energy Eq in the unit time 
A t, the total Fe Ka emission is given by 



L{FeKa) = nL'{FeKa) 



n 



So 1 E^<'Y, 



FeKa^FeKa 



At hp 



K 



gas 



Rn 



(2) 



where and Yp^^^ are the absorption opacity and the FeKa yield of z-times ionised iron, 
^^FeKa ^hc energy of the Ka line of i-times ionised iron (~6.4keV) and the summation 
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is over all abundant ionisation stages, k^J^'^ is the absorption opacity due to all other abun- 
dant species and -R^. is the branching ratio between Ka and K[3 fluorescence (0.882:0.118, 
Bambynek et al.l Il972l). We use a value of 0.34 for the fluorescence yields of neutral and 
once-ionised iron (iBambynek et al.l Il972l : iKrausd 1 19791 : iKallman et al.l 120041 ) . 



The n FeKa energy packets are emitted in random directions as the fluorescent emis- 
sion process can be assumed to be isotropic. To ensure conservation, energy per unit time 
{eo/At — L{Fe Ka)) must also be immediately re-emitted from the event location, at a fre- 
quency determined by s ampling the local emis sivity of the medium and transferred through 
the grid as described by lErcolano et al.l (120031 ). 



The fates of the newly emitted FeKa packets are determined by the absorption and 
Compton opacities encountered along their diffusion paths. Fe Ka packets that are absorbed 
will be transformed into diffuse field packets or emission line packets at another energy and 
will not contribute to the emerging spectrum at 6.4 keV. FeKa packets that are Compton 
scattered will experience a change of direction and a shift in frequency (by approximately 
(/iz/)^/meC^), their new frequency and direction being determined using the Klein- Nishina 
formulae. These packets may therefore finally emerge and contribute to the low energy 
shoulder of the Fe Ka feature, or may undergo further scatterings and/or absorptions. 

The emergent integrated and the direction-dependent spectral energy distributions (SEDs) 
are finally calculated using all escaped energy packets, including contributions from Compton 
reflected packets and fluorescence hne packets. Our method also allows us to easily separate 
the various components of the emergent SED. 

Since the study presented here is similar to the earlier B79 work, we summarise in brief 
the differences between these studies: 



1. B79 was aimed at understanding Fe Ka fluorescence in solar flares; the possible utility 
of fluorescence for estimating the photospheric Fe abundance (which was of course 
more uncertain that at present) was emphasised. The calculations presented here 
extend the parameter space so as to cover conditions found in large stellar flares and 
stars completely covered in X-ray emitting active regions; our emphasis is on the use of 
fluorescence as a diagnostic of coronal and flare geometry on distant unresolved stars. 

2. Input parameters used here are more up-to-date in several respects. Input spectra have 
been computed using isothermal coronal plasma models whereas B79 used an analytical 
approximation to the bremsstrahlung continuum. B79 also adopted a simple form for 
absorption cross-sections of Fe and the photospheric gas. Here, photoabsorption cross- 
sections are computed self-consistently according to the chemical co mposition and ioni- 
sation state of the gas using analytical fits to Opacity Project data fIVerner fc Yakovlev 
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19951 ). These abundance- weighted cross-sections are about 40% higher than that used 
by B79 for the gas with Fe excluded; the B79 cross-section for neut ral Fe K-shell ab- 
sorption is about 20% larger than that of IVerner fc YakovlevI (Il995l ) at threshold, but 
very similar at energies above 8 keV. As we shall see below, these differences have little 
influence on the resulting fluorescence efficiencies, largely because the Fe cross-section 
dominates and differences here are small. 



3. B79 assumed that once a Fe Ka line photon is Compton scattered it is "lost" due to 
the shift in energy. Here, we follow all photons until they emerge from the photosphere 
and therefore also resolve the Compton shoulders comprising scattered line photons. 



2.3. Benchmark Calculations of Iron Ka in solar flares 

In order to establish the robustness of our new code, we first examine the Fe Ka problem 
for point-like solar flares using the full 3D geometry of our models and fully self-consistent 
radiative transfer, treating photoabsorption contributions from all abundant atoms and ions, 
and (Compton) scattering by electrons. As a benchmark, we attempt to reproduce the earlier 
results of B79. 

The photon density spectrum of the incident X-rays, L{e), is taken to follow the same 
bremsstrahlung power-law with energy, e, as defined in B79: 

L{e) = A e-^ exp{-e/kT) photon s'^keV"^ (3) 

We consider four X-ray temperatures (0.5, 1.0, 3.0, 5.0 keV) and two flare heights, h = 
and h = 0.1 Rq, which cover the full parameter space investigated by B79. In the remainder 
of this paper we will quote temperatures in degrees K, however in this section these are 
given in keV, which are the units used by B79. For the sake of comparison, we used the Fe 
abundance adopted by B79 of npe/nn = 4 x 10^^ by number; for other elements we adopted 
the mixture of GS. In Table [1] we compare our results for the iron Ka fluorescence efficiency, 
F, to those of B79 (values in brackets). F is defined as the ratio of the total luminosity of 
iron Ka photons emitted from the photosphere and A^7.ii/2, where Nj u = f!^^-^L{e)dE, 
with L{e) given by Equation [31 



In this formulation the flux of Fe Ka photons received at Earth is equal to 

r/(^)iV7.ii 



Ka 



(4) 



47rD2 ' 

where D is the distance from the source and 6 the heliocentric angle subtended by the flare 
and observer. 
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The agreement with the results of B79 is generally very good, and significantly better 
at higher X-ray temperatures; discrepancies amount to only 18% for X-ray temperatures of 
0.5 keV. We attribute the bulk of the differences in fluorescent efficiencies to the use by B79 
of slightly different absorption cross- sect ions and an analytical approximation for the gas 
opacities, while we calculate these self-consistently according to the elemental abundances 
considered in our models. The slight temperature dependence in the differences is most likely 
due to the different slopes of the B79 functional form and our calculations, which take into 
account all elements with Z<30. 

Our 3D models allow us to obtain the spectral energy distribution emerging from arbi- 
trary viewing angles. In Figure [3] we show the variation of f{6) for ten heliocentric angles 
from 5° to 95° for two flare heights {h = and h = O.IRq). These plots are directly compa- 
rable with the top-left and bottom-right panels of Figure 3 of B79 and show good agreement 
with his results. 

Finally, in Figure H] we show the variation of F as a function of iron abundance. Again, 
our results are in good agreement with those of B79 (c/his Figure 4) and show that the 
relation of the FeKa fluorescence efficiency on the photospheric iron abundance is weaker 
than a proportional relation, as expected. We discuss the Fe dependence of the Fe Ka flux 
in more detail in §3.41 below. 



3. New Fe Ka Diagnostic Calculations 



Having in the previous section established the robustness of our methods, we present the 
results from a new grid of model calculations designed to provide useful diagnostics for FeKa 
detections for a wider range of stellar environments. In particular, stellar coronae cover a 
much wider range of temperatures and Fe abundance than the solar case and investigation 
of fluorescence over a wider range of parameters will be necessary for correct interpretation 
of fluorescent lines from stars. 

We investigate two geometrical configurations: (i) a single flare and (ii) an active corona 
that can be approximated by a spherical shell illuminating a photosphere from particular 
scale height h. While coronae are of course not expected to conform to a shell-like geometry, 
this case should approximate quite well the case in which a star is quasi-uniformly covered 
in active regions. Coronal fluorescing spectra were adopted from a grid of isothe rmal models 



computed using emissivities from the CHIANTI compilation of atomic data (jLandi et al. 



20061 . and references therein), together with ion pop ulations from iMazzotta et al.l (Il998l ). as 



implemented in the PINTofALE IDL software suit flKashvap fc Drake! l20o"o[ l We adopted 
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the chemical composition of GS for all calculations except where noted. While there will 
be some contribution from transitions in He-like and H-like Ni and the higher Lyman series 
of He-like and H-like Fe to the ionising photon flux A^7.ii, this is small compared with the 
integrated continuum contribution; consequently the metallicity sensitivity for the ionising 
spectrum is generally negligible for realistic ranges of metallicity. 

3.1. A Single Flare Illuminating a Photosphere 

In the case of a single flare illuminating a photosphere from a given height, the heliocen- 
tric angle, 9, has a dramatic effect on the FeKa flux detected. The intrinsic FeKa efficiency, 
r, is, of course, independent of 6*, being a function of the fluorescing spectrum, the flare 
height and the relative iron abundance. Using GS solar abundances for all elements with 
Z < 30 we show in Table [2] the FeKa efficiency, F, predicted for X-ray temperatures in the 
range T = 6.4-8.0 keV and flare height in stellar radii in the range h = O-Si?,,. 

In order to investigate the Fe Ka line strength as a function of heliocentric angle, we 
calculated f{6) for a number of angles at each flare height. For convenience we fit our results 
using 4th or 5th order polynomials, such that f{6°) = Yl^=o Q^i " The fit coefficients are 
listed in Tableland the resulting f{9) functions are plotted in Figure O Our fits are valid 
between heliocentric angles 0° and 100°; the maximum errors between our data and the fits 
at each flare height are also listed as percentages in Tableland are less than 2.4%. 

These functions illuminate interesting characteristics of the fluorescence problem for 
different incidence and exit angles. Firstly, we point out that for flnite flare heights h, f{6) 
is non-zero for angles 6 > 90°. This is simply due to the surface illumination from flares 
beyond the limb extending onto the visible hemisphere: the angle at which f{6) reaches is 
given geometrically by Oq = 90° + arccos{R^/ [Ri, + h)). We note, however, that we did not 
explore the very low-response tail of f(9) for heliocentric angles greater than 100°, as the 
Fe Ka flux, although flnite for some flare heights, is too small to be detectable even for very 
strong flares. Furthermore, f{9) falls below the statistical variance of our models at these 
large angles. 

The general decrease of f{9) with increasing 9 present in all curves is largely a result of 
the l/cos{9) increase in path length required for Ka photons to escape. The most conspic- 
uous characteristic of f{9) behaviour with different heights is the steepening of the decline 
from smaller to larger angles with increasing flare height. This arises because of the range of 
incidence angles of flare photons on the photosphere. As h ^ oo, incidence angles relative to 
the photospheric normal become smaller and the average penetration depth before interac- 
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tion by either absorption or scattering events is larger. Consequently, path lengths increase 
and Fe Ka escape probabilities decrease. Folded in with this is the decreasing fraction of 
photons incident on the visible hemisphere of the star with increasing 6. 



3.2. A Coronal Shell Illuminating a Photosphere 

The calculations for f{9) and F can be readily adapted from the point source case to 
one in which the corona is approximated by a spherical shell of emission, as might be the 
case for stars significantly more active than the Sun. This can be seen by taking the shell 
case to be a large ensemble of point sources distributed uniformly about the star for which 
the observed fluorescent flux is then simply the solid angle average of the point source case, 

^^-^ ' ^mdn (5) 



4nD^ 47r 

Since, following from the definition of F and our Equation 4, J f{6) dQ = 2it, the flux is 

By tailoring the solid angle integral in Equation [5] the fluorescent flux for an arbitrary 
coverage pattern of active regions viewed from any angle can be derived from the tabulated 
values of f{9) and F. 



3.3. Sensitivity to Ionising X-ray Temperature 

The fluorescence efficiency, F, is illustrated as a function of the isothermal plasma tem- 
perature adopted for the ionising coronal spectrum in Figure El We find F varies very slowly 
with temperature, decreasing by only a factor of ~ 2 over a factor of 30 or more in tem- 
perature. This gradual decrease is predominantly a result of the competition between Fe 
inner-shell absorption and Compton scattering. As temperatures increase, incident photons 
are present at higher and higher energies. Since the Fe cross-section declines with increasing 
photon energy while the Compton cross-section remains constant, higher energy photons are 
more liable to be scattered out of the photosphere before being destroyed by photoabsorption. 



3.4. Sensitivity to Photospheric Fe Abundance and Metallicity 



In the optically-thick fiuorescence case such as characterizes the photospheric fiuores- 
cence problem, inner-shell photoabsorption by Fe atoms and ions competes with the back- 
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ground opacity from other elements. Observed Fe Ka line fluxes and equivalent widths are 
then expected to vary predominantly according to the relative Fe abundance rather than just 
the photospheric metallicity. Note that the optically-thick situation differs to the optically- 
thin formalism discussed by Liedahll ( 19991) arid referred to in the discussion of stellar and 
protostellar Fe Ka lines (e.g. lFavata et al.ll2005l : lTsujimoto et al.ll2005l : lOsten et al.ll2007l ): in 
the optically-thin case the fluorescent Fe line strength is primarily a function of the absorbing 
Fe column. 



In extension to the benchmark calculations discussed in §2.31 we have computed the 
sensitivity of Fxa to the Fe abundance for a model fluorescing spectrum with a temperature of 
10^-2 K and Fe abundance in the range -0.8 <[Fe/H]< +0.3 (0.5x10"^ < npjn-a < 6x10"^). 
This temperature was chosen to be representative of a typical very active corona or moderate 
stellar flare. Trends with [Fe/H] are expected to be only weakly sensitive to this adopted 
temperature, the fluorescence efficiency, F, varying only very slowly with T (Figure [6]). The 
results are illustrated in Figure [7] where we show F normalised to its value at the GS solar 
Fe abundance in comparison to the proportional relation, analogous to Figure H] discussed in 
gSJland Figure 4 of B 79. 



The fluores cence line s t rength vs Fe abundance was investigated in the context of ac- 
cretion disks by lMatt et al.l (119971 ). and our results are comparable to Figures 1 and 2 from 
that work for the range of Fe abundance in common. As noted above, B79 showed that the 
relationship between the relative fluorescence efficiency F and the Fe abundances is weaker 
than a proportional relation. In fact, the relation increasingly departs from proportionality 
with increasing Fe abundance. This behaviour arises because of the equipartition of photoab- 
sorption between the different species present in the atmosphere. At low Fe abundance when 
the background cross-section due to elements other than Fe dominate, F^a varies according 
to npf.apeK/o'tot, where ap^K and atot are the Fe n = 1 shell and total plasma absorption 
cross-sections, respectively. At very large Fe abundance, the photoionization cross-section 
of Fe will begin to dominate atot] crpeK/o'tot then tends to asymptote to the constant ratio 
0"FeK/(c"FeK + ctfbL.m), where (TFeL.M is the total Fe photoabsorption cross-section from levels 
n > 1. We do not consider such large Fe abundances here; this asymptotic behaviour is nicely 
illustrated by lMatt et al.l (119971 ). The relationship is further complicated to some extent by 
the role of Compton scattering, but this becomes less relevant for higher Fe abundances for 
which photoabsorption dominates Compton scattering out to higher energies. 

Also illustrated in Figure [7] is the relative fluorescence efficiency for different pho- 
tospheric metallicity (ie where all metal abundances are scaled together) for the range 
—1.0 <[M/H]< +0.3 relative to the GS solar mixture. As expected, the fluorescent effi- 
ciency is not so sensitive to [M/H] because of the tendency of Fxa to follow the relations 
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above between the different opacity sources: adjusting the metalhcity as a whole results 
in a smaller change in npeCTFeK/ctot since atot is by definition in lockstep with metallicity. 
As metallicity decreases, Compton scattering and the background H and He cross-sections 
become more important, the combined effect of which is to weaken F^a in comparison to 
the higher metallicity case. For very high metallicities, these opacity sources are negligible 
and the fluorescence efficiency asymptotes to the value given by npeO'FeK/o'tot- 



4. Discussion and Applications 

The main scientific motivation for this work is to provide the foundation to use Fe 
fluorescence as a quantitative diagnostic of coronal and flare geometry. There now exists a 
handful of detections of fluorescent emission from stars. Sensitivity is currently limited to a 
large extent by the low spectral resolution of available instruments and progress is expected 
to accelerate dramatically with the future availability of X-ray calorimeters. 

Since modern X-ray spectral analyses based on low-resolution CCD pulse-height spectra 
tend to express line strengths in terms of the line equivalent width, we have computed this 
quantity for the case of ^ = and the ranges of heights and X-ray temperatures investigated 
in §ni The equivalent width in this context refers to the fluorescent, processed line seen on 
top of the continuum of the ionising coronal spectrum. While the 6 = case gives the most 
optimistic line strength, we note that f{6) is quite slowly varying for angles 6 < 45° for 
the flare heights for which significant FeKa might be observed. The equivalent widths are 
illustrated in Figure [8] and listed in Table HI 



4.1. Fluorescence from Pre-Main Sequence Stars 

The observability of the cold Fe Ka line is of course strongly dependent on the quality 
of the X-ray spectrum obtained. The most extensive study of PMS Fe fluorescence to date is 



that based on Chandra observations of the Orion Nebula Cluster by lTsujimoto et al.l (120051 ). 
This study detected significant 6.4 keV excesses attributable to Fe fluorescence for 7 out of 
127 sources found to have significant counts in the 6-9 keV band. Equivalent widths were 
in the range 110-270 eV at plasma temperatures of ~ 3-10 keV. There is clearly a strong 
selection effect here and these fiuorescent line strengths likely represent the upper end of the 
distribution. 



Our calculations for a flare at scale height h = Oi?^ are also appropriate for a flare occur- 
ring above an inflnite plane, such as might approximate a disk-encircled PMS star. As in the 
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photospheric case, the fluorescence problem can be treated orthogonally from the ionisation 
structure of the disk, which is not greatly altered from its very largely neutral overall state 
by X-rays from a typical flare. Any small degree of X-ray photoionisation will also not affect 
Fe Ka line strengths because fluorescence yields are essentially invariant for lower Fe ions. 
Our calculations indicate that attaining equivalent widths much in excess of 100 eV is not 
straightforward for such a simple geometry for the plasma temperatures observe d in the flu 



oresci ng X- ray spectra. This finding is in agreement with earlier calculations by iMatt et al. 
( I1991I ) and lGeorge fc FabianI (Il99ll ). who find equivalent widths of ~ 150 eV for a flat disk 
illuminated by X-rays with power-law photon spectral energy distributions. 

There are at least four ways in which equivalent widths might be elevated above the 
values we find: (1) super-solar Fe abundance in the disk material, possibly arising as a result 
of an elevated dust-to-gas ratio; (2) disk flaring, resulting in a solid angle coverage > 2n; (3) 
line-of-sight obscuration of the central flaring source (but not fluorescent line photons) by 
optically thick structures such as the star itself (ie the flare occurring on the far hemisphere); 
and (4) fluorescence contributions from ionisation by non-thermal electrons. 

By an alogy with the solar case, in which exci tation by non-thermal electrons is usually 
negligible (jParmar et al.l Il984l : lEmslie et al.l 119861 ) and i s much more difficult on ene rgetic 
grounds, we consider (4) the least plausible of these. iBallantyne fc Fabian! (120031 ) have 
also shown in the accretion disk context that Fe K production by non-thermal electron 
bombardment requires 2-4 orders of magnitude greater energy dissipation in the electron 
beam than is required for an X-ray photoionization source. 

Disk fiaring can give rise to increased line strengths by factors < 2 simply from con- 
sideration of the increased solid angle coverage possible compared with an infinite fiat disk. 
It is also difficult to envisage enhanced Fe abundances in the disk being able to elevate 
line strengths by more than a factor of a few. Of some interest, then, is the observation 
of an Fe Ka line equivalent width of ~ 1400 mA during t he rise phase of a fiare on the 
PMS Orion nebula star forming region object V 1486 Ori by lCzesla fc SchmittI (120071 ). Such 
an enhancement over an infinite disk value of ~ 150 mA, even with a large degree of disk 
flaring, would still require extrem e enhancenients o f the disk Fe abundance by an order of 
magnitude or more (e.g. §3.41 and iMatt et al.lll997l ) were the line due to photoionisation by 
the directly observed continuum. We point out, however, that fluorescent line photons from 
a PMS disk can still be observed when the X-ray flaring source is located behind the star 
and obscured from the line-of-sight. In the case of the V 1486 Ori flare, the large observed 
Fe Ka equivalent width is simply and plausibly associated with a partially obscured flare 
whose rise phase was not fully observed directly owing to line-of-sight obscuration by the star 
itself. Such obscured flares will inevitably be the cause of some fraction of observed Fe Ka 
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lines from PMS stellar disks. This explanation is also more consistent than one relying on 
preferential disk geometry and Fe abundance with the non-detection of Fe Ka from a second 
less extreme flare whose impulsive phase instead appears to have been quite visible. 



4.2. Fluorescence from stellar photospheres 



The strong flare in II Peg observed by Swift and analysed by lOsten et al.l (120071 ) presents 
another interesting case. Photospheric Fe Ka was clearly detected throughout the event. 
Equivalent widths for different times in the flare ranged from 18 to 61 eV, with uncertainties 
in the 20-45% range. The authors favoured a coUisional excitation mechanism for the line, 
arguing that fluorescence would be unlikely to produce an observable feature. This assess- 
ment employed a simple analytical formula applicable to optically-thin cases in which only a 
minor fraction of the incident X-ray flu x is subject to p hotoa bsorption or scattering (see e.g. 
Liedahll EqQ^ : iKrolik fc KallmanI Il987h . lOsten et all J2007h correctly noted that the path 
length required to obtain the observed equivalent widths under such an approximation was 
similar to the r = 1 Compton scattering depth, but discounted fluorescence as a possibility 
on these grounds. Other than the inapplicability of the optically-thin formula for the photo- 
spheric fluorescence case, one reason such an argument is overly pessimistic is that incidence 
angles on the photosphere range from ~ 0-90° for small scale heights and path lengths for 
escape are less than penetration depths by the factor of the inverse cosine of these angles. 
We defer a more detailed treatment of the event to future work, but note here that equivalent 
widths of 50 eV are achieved for flare heights up to ~ 0.2R^ or so for the 10* K model in our 



grid, a temperature similar to the average of the values found for the flare by lOsten et al. 



(120071 ). While collisional ionisation cannot be ruled out observationally as the source of the 
observed Fe fluorescence, it is not a requirement. 



5. Summary 

We have investigated the production of Fe Ka fluorescence lines by irradiation from 
coronal and flare emission using a 3-dimensional Monte Carlo radiative transfer approach 
including Compton redistribution. The results are presented in the form of convenient tables 
describing the fluorescence efficiency as a function of the flare height, the temperature of the 
ionising X-ray spectrum, and the viewing angle. We have also illustrated how the fluorescent 
efficiencies scale with the photospheric metallicity and the Fe abundance. The results should 
be of use for interpreting observations of Fe Ka lines seen from stars. 
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Our computed Fe Ka equivalent widths for irradiation from a height of zero above the 
photosphere correspond to the case of flares above a disk of infinite extent and are relevant to 
observations of fluorescence from PMS stars. Observed equivalent widths tend to be slightly 
larger than our computed ones for a plasma of solar composition and are an indication of 
flaring disk geometries or super-solar Fe abundances. For one extreme case recently observed 
on V 1486 Ori, we propose that the very large equivalent width observed (~ 1400 mA) arose 
from X-ray fluorescence by a flare partially obscured from the hne-of-sight by the stellar 
limb. 

The FeKa equivalent width reported for a large flare on 11 Peg is consistent with our 
computed values for a flare scale height of a few tenths of a stellar radius, and a coUisional 
excitation mechanism is not a requirement. 
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Tables 



Table 1: Fluorescence efficiencies, F [%], for B79 benchmark models. B79 values are given 
in parentheses. The iron abundance ripe = 4 -lO^^nn was used. 







X-ray temperature (keV) 




Height (Re) 


0.5 


1.0 3.0 


5.0 


0. 


4.06 (3.44) 


3.70 (3.26) 2.93 (2.74) 


2.47 (2.44) 


0.1 


1.99 (1.78) 


1.80 (1.65) 1.39 (1.37) 


1.16 (1.22) 
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Table 2: Fluorescence efSciencies, T [%], for a range of X-ray temperatures, Tx, and flare 
heights. iGrevesse &: Sauvall (119981 ) (GS) solar abundances were adopted for all elements. 





X-ray temperature, 


Log(Tx; 


[K] 


Height (R^) 


6.4 


6.8 


7.2 


7.6 


8.0 


0.00 


3.95 


3.38 


2.87 


2.16 


1.69 


0.05 


2.70 


2.27 


1.84 


1.58 


1.22 


0.10 


2.12 


1.78 


1.45 


1.23 


0.96 


0.50 


0.76 


0.63 


0.52 


0.44 


0.34 


1.00 


0.44 


0.30 


0.25 


0.20 


0.16 


5.00 


0.047 


0.042 


0.032 


0.026 


0.021 



Table 3: Polynomial fit coefficients to heliocentric angle functions. The maximum error 
between our fit and the data is given in the last common as a percentage value. 



Height (R,) 




ai 


a2 


«3 


^4 


"5 






0.00 


1.45 


-1.53E-3 


2.92E-4 


-1.73E-5 


2.88E-7 


-1.68E-9 


0.6 


0.05 


1.49 


1.26E-3 


-3.97E-5 


-1.40E-6 


-2.54E-9 


0. 


1.6 


0.10 


1.56 


5.78E-3 


-6.37E-4 


2.02E-5 


-3.07E-7 


1.45E-9 


1.2 


0.50 


1.75 


1.90E-3 


-2.17E-5 


-6.37E-6 


5.09E-8 


6.63E-11 


2.3 


1.00 


1.88 


5.36E-3 


-3.99E-4 


3.46E-7 


l.OlE-8 


0. 


0.8 




5.00 


2.25 


9.00E-3 


-1.75E-3 


4.19E-5 


-4.69E-7 


1.94E-9 


2.3 



Table 4: FeKa line e quivalent widths (EW) at 6 = for a range of X-ray temperatures, 

Tx, and flare heights. iGrevesse fc Sauvall (119981 ) (GS) solar abundances were adopted for all 

elements. 

X-ray temperature, Log(Tx) [K] 



Height (R,) 


6.4 


6.8 


7.2 


7.6 


8.0 


0.00 


0.54 


6.73 


17.73 


59.45 


132.57 


0.05 


0.39 


4.79 


14.71 


49.49 


98.66 


0.10 


0.33 


3.88 


12.73 


39.95 


78.23 


0.50 


0.14 


1.63 


6.02 


18.47 


34.87 


1.00 


0.076 


0.86 


3.28 


9.92 


17.52 


5.00 


0.0054 


0.12 


0.46 


1.23 


2.22 
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Fig. 1. — The photoabsorption cross-section in the energy range of interest for a neutral gas (sohd curve) 
and once-ionised (dashed curve) plasma of solar composition. The labels indicate the ionisation edges for 
most abundant elements. 
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Fig. 2. — The locus corresponding to optical depth of unity for X-rays incident on the solar photospheric 
VALC model plotted as a function of incident photon energy and height (thin black curve). Toward lower 
energies absorption occurs predominantly in the lower chromosphere; this region is shifted downwards with 
increasing energy, reaching the upper photosphere at the Fe K photoionization threshold energy. The model 
temperature structure as a function of height is illustrated by the thick grey curve. 
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Fig. 3. — The dependence of the FeKa fluorescence flux on the heliocentric angle, 
parametrized by the f{6) factor, for two heights. These results are directly comparable 
to the top left and bottom right panels of Figure 3 of B79. 
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Fig. 4. — Dependency of FeKa efficiency on the abundance of iron for the B79 benchmark 
case. The proportional relation is indicated by the dashed fine. This figure is directly 
comparable to Figure 4 of B79. 
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Fig. 5. — Fits to the f{6) factor relating the FeKa fluorescence flux to the angle of obser- 
vation obtained for six flare heights 
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log T [K] 

Fig. 6. — The fluorescence efliciency, F, as a function of the ionising X-ray spectrum tem- 
perature for different flare heights, computed for the GS solar abundance mixture. 
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Fig. 7. — The fluorescence efliciency, F, as a function of logarithmic photospheric metalhcity 
[M/H] and Fe abundance [Fe/H], normahsed to its value for the solar mixture of GS. The 
efficiencies as a function of [Fe/H] were computed for GS values for the abundances of all other 
elements. The dashed curve represents the proportionality relation expected from scaling 
the efficiency linearly with Fe abundance from the lower Fe abundance point. Calculations 
where performed for an ionising X-ray spectrum with temperature of 10^'^ K. 
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Fig. 8. — The Fe Ka fluorescent line equivalent width at = as a function of X-ray 
temperature and flare height, computed for the GS solar abundance mix. 



